Athermal Shear- Transformation-Zone Theory of Amorphous Plastic Deformation II: 
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In the preceding paper, we developed an athermal shear-transformation-zone (STZ) theory of 
amorphous plasticity. Here we use this theory in an analysis of numerical simulations of plasticity in 
amorphous silicon by Demkowicz and Argon (DA). In addition to bulk mechanical properties, those 
authors observed internal features of their deforming system that challenge our theory in important 
ways. We propose a quasithermodynamic interpretation of their observations in which the effective 
disorder temperature, generated by mechanical deformation well below the glass temperature, gov- 
erns the behavior of other state variables that fall in and out of equilibrium with it. Our analysis 
points to a limitation of either the step-strain procedure used by DA in their simulations, or the STZ 
theory in its ability to describe rapid transients in stress-strain curves, or perhaps to both. Once 
we allow for this limitation, we are able to bring our theoretical predictions into accurate agreement 
with the simulations. 
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I. INTRODUCTION 

In the preceding paper we presented an athermal 
STZ theory of plastic deformation in materials where 
thermal activation of irreversible molecular rearrange- 
ments is negligible or nonexistent. Here we use that the- 
ory to interpret an extensive body of computational data 
published recently by Demkowicz and Argon 0, d, 0, S] , 
hereafter referred to occasionally as "DA." Those authors 
simulated plastic deformation in amorphous silicon using 
a system of 4096 atoms interacting via a Stillinger- Weber 
potential Q in a cubic cell with periodic boundary con- 
ditions. They subjected this system to pure shear under 
both constant-volume and constant-zero-pressure, plane- 
strain, conditions. Their reports of these simulations are 
remarkably complete and detailed. They provide valu- 
able and challenging information about: 

• the relation between stress-strain response and 
sample preparation; 

• the theoretical description of nonequilibrium be- 
havior in systems subject to steady-state and tran- 
sient mechanical deformation; 

• the nature of the glass transition in simulated amor- 
phous silicon; 

• the strengths and limitations of numerical simula- 
tion techniques. 

We address each of these topics in the body of this report. 

Demkowicz and Argon used two different procedures 
for simulating shear deformation. In their potential en- 
ergy minimization method (PEM), each step in the pro- 
cess consisted of a small, affine, shear displacement of 
all the atomic positions, followed by a minimization of 
the potential energy during which the atoms relax to 
their nearby positions of mechanical equilibrium. Sup- 
posedly, PEM simulations correspond to the limit of zero 
strain rate at zero temperature, but that interpretation is 



problematic. In their molecular dynamics (MD) method, 
DA used a different step-strain procedure in which each 
small shear increment was followed by an MD relaxation 
at temperature T = 300 iC, with an average strain rate 
of order lO^sec"^. In both procedures, an incremental 
shear was imposed only after the system was judged to 
have reached a stable, stationary state following the pre- 
ceding step. We argue below that there are important 
uncertainties associated with both of these step-strain 
simulation methods. 

The single most important feature of the DA simula- 
tions is that, in addition to measuring the shear stress 
(and keeping track of pressure and/or volume changes) 
during deformation, DA also observed local atomic corre- 
lations within their system. Here they were taking advan- 
tage of their numerical method to see inside their system 
in a way that is seldom possible in laboratory experi- 
ments using real materials. They found that the environ- 
ments of some atoms were solidlike and others liquidlike, 
and that the liquidlike regions seemed to be, as they say, 
the "plasticity carriers." Before any mechanical deforma- 
tion, their fraction of liquidlike regions (f> was small when 
the system was annealed or cooled slowly, and was ap- 
proximately 0.5 when the system was quenched rapidly. 
Then, during constant-zero-pressure deformations, ((> ap- 
proached a value slightly less than 0.5, independent of 
its initial value. Thus, 4> behaved in a manner qualita- 
tively similar to the dimensionless effective temperature 
X or, equivalently, the density of STZ's A oc exp(— 1/x) 
described in [l]. The relation between these quantities is 
one of the main topics to be addressed in this paper. 

Amorphous silicon, like water, expands as it solidi- 
fies. Moreover, its properties are highly sensitive to small 
changes in density. A slowly quenched, more nearly equi- 
librated system is less dense than one that is rapidly 
quenched, because the former contains a smaller popula- 
tion of denser, liquidlikc regions than the latter. When 
strained, the more nearly equilibrated system initially 
responds elastically, and takes longer than a rapidly 
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quenched system to generate enough plasticity carriers 
to enable plastic flow. Accordingly, the equilibrated sys- 
tem exhibits a more pronounced stress peak of the kind 
illustrated in Figure 1 (top curve, upper panel) of [l|. As 
the liquidlikc fraction (f) increases, the system contracts 
if held at constant pressure; or else, if the system is held 
at constant volume, the pressure decreases and may even 
become negative. It may be surprising to theorists but 
nevertheless is true that, in this material, the STZ's must 
be associated with decreases rather than increases in free 
volume. 

For theoretical purposes, it will be simplest to look only 
at constant-pressure simulations, because the changes 
in pressure that occur at constant volume require addi- 
tional, hard-to-control approximations for behaviors that 
are not intrinsic to the STZ hypotheses. System-specific 
parameters such as the yield stress Sy are likely to be 
more sensitive to changes in pressure than in volume. 
For example, the Mohr-Coulomb effect implies that Sy 
increases with pressure. Demkowicz and Argon point out 
that direct evidence for pressure dependence can be seen 
in the graphs of the liquid fraction f6 as a function of 
strain 7 in Figs. 10(c) and 15(c) of Q. In the constant- 
zero-pressure simulations shown in their Fig. 15, <j) equi- 
librates to approximately the same value after deforma- 
tion for all four different initial conditions. That does 
not happen in the constant-volume simulations in Fig. 
10, implying that the internal states of these deformed 
systems differ from one another in nontrivial ways. 

In addition to the questions pertaining to step-strain 
procedures, to be discussed in detail in Sec. IIVI there 
are other limitations of the DA simulations that must be 
taken into account. The DA simulation system is too 
small to make it likely that more than one STZ-triggered 
event is taking place within a characteristic plastic relax- 
ation time; and data is reported only for single numerical 
experiments performed on individual samples rather than 
averages over multiple experiments. Thus, fluctuations 
are large, and the results must be sensitive to statisti- 
cal variations in initial conditions. In short, we must be 
careful in our interpretations. 

On the positive side, Demkowicz and Argon's remark- 
able combination of multiple simulation methods and 
multiple observations allows us to use the athermal STZ 
theory to construct what we believe is an internally self- 
consistent interpretation of their data. Going beyond 
their stress-strain curves, and focusing on the behaviors 
of the liquidlike fraction (j) and the density p, we develop a 
quasithermodynamic picture in which the effective disor- 
der temperature x plays a dominant role. In this picture, 
the STZ's, i.e. the active flow defects, are rare sites, out 
in the wings of the disorder distribution, that are more 
susceptible than their neighbors to stress-induced shear 
transformations. They almost certainly lie in the liquid- 
like regions of the system, but only a very small fraction 
of the liquidlike sites are STZ's. The energy dissipated 
in the STZ-initiated transitions generates the effective 
temperature Xi a systemwide intensive quantity that, in 



turn, determines extensive quantities such as 4> and p. 

Our central hypothesis is that there exist quasithermo- 
dynamic equations of state relating steady-state values 
of (j) and p to X (and also, in principle, to the pressure 
and shear stress). Once we have found these equations 
of state, we extend the quasithermodynamic model to 
describe the way in which 4> and p fall in and out of 
equilibrium with x during transient responses to exter- 
nal driving forces. In this way, we arrive at a quantitative 
interpretation of the DA simulations. 

We begin in Sec. [TT]with a brief summary of the ather- 
mal STZ theory developed in [ij . Section IIIII contains 
a preliminary analysis of the DA stress-strain data. In 
Sec. IIV[ we explain why we cannot accept the initially 
good agreement between the STZ theory and the simula- 
tion data. These arguments motivate the quasithermody- 
namic hypothesis, introduced in Sec. [V] We extend the 
quasithermodynamic ideas to nonequilibrium situations 
in Sec. IVII Finally, in Sec. IVIIl we conclude with some 
speculations concerning the validity of the quasithermo- 
dynamic picture and its relation to the STZ theory in 
general. 

II. STZ SUMMARY 

In order to make this paper reasonably self-contained, 
we start by restating Eqs. (3.12) - (3.16) in [l|. The first 
of these equations is an expression for the total rate of 
deformation tensor = 7/2 as the sum of elastic and 
plastic parts: 

i5*°* = -i^+i?f'(J,m,A), (2.1) 
2 p 

where s and p, are the deviatoric stress and the shear 
modulus measured in units of the yield stress Sy. The 
internal state variables m and A are, respectively, the 
orientational bias and the scaled density of STZ's. We 
consider only the case of pure, plane-strain shear in which 
the material is strained at a fixed rate 7, and the stress is 
measured as a function of the strain 7. To describe such 
experiments, we write Eq. (|2.ip in the form 

ds , ( 2eo A ,^ A 

— =^ l--i^g(s,m) , (2.2) 

where eo is a number of order unity (the atomic density 
in atomic units multiplied by the incremental strain as- 
sociated with an STZ transition), tq is the characteristic 
time scale for STZ transitions, and 

g(s,m)=C(s)(^-m). (2.3) 

The function C(s) describes the stress dependence of the 
STZ transition rate; it is proportional to |s — 1| for |s| ^ 1 
and vanishes smoothly near s = 0. As defined by Eq.(3.6) 
in [1], the smoothness of C(s) near s = is determined 
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by a parameter which we take to be unity throughout 
the following. 

Similarly, we restate Eqs. (3.9), (3.10) and (3.11) from 

^ = J-Sq(S^m)(e-^^^ ~A); (2.5) 
d-y jtq 

and 

^= Asq{s,m) (xoc - x)- (2-6) 

c(7 Co 7 To 

The constant cq is a configurational specific heat per 
atom in units ks , which must be of order unity. 

We assume that the stress and strain tensors in these 
equations remain diagonal in two dimensional pure shear 
and plane strain, and that any correction in the third 
dimension is negligible in comparison to other uncertain- 
ties in this analysis. A basic premise of the STZ theory 
is that the zones are rare and do not interact with each 
other; thus we assume from the beginning that the quan- 
tity eo exp(— 1/x) is small of order 10"'^ (in fact, very 
much less), so that the equations for m(7) and A(7) are 
stiff compared to those for 5(7) and xil)- Then we can 
safely replace Eqs. (|2.4p and (|2.5p by their stationary 
solutions: 

;f!!|;; <"> 

and 

A = e~i/^. (2.8) 
Eqs. dMl) and dM]) become: 

^-Mfl-^e-Vx,(,)), (2.9) 

and 

^ = lfle-i/'^Sg(S)(Xoo-x), (2.10) 
d-y Co go 

where 

(7(l)-C(S)(^-mo(S)), (2.11) 

and qo = 7 To. 

III. FIRST FITS TO THE DA STRESS-STRAIN 
CURVES 

The first step in our STZ analysis of the DA data is to 
use Eqs. (|2.9p - (|2.1ip to fit the stress-strain curves in Fig. 
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FIG. 1: Theoretical stress-strain curves (solid lines) compared 
to the DA numerical simulation data [31 . The parameters 
used are: eo = 1, s f — 1.06GPa, /i = 46 GPa, Sy = 0.35 GPa, 
Co = 0.18 and Xoo = 0.065. The initial effective temperatures 
are: xo= 0.0559, 0.0580, 0.0612, 0.0680 (from top to bottom). 
The curves are shifted by 0.5 GPa for clarity. 

15(a) of [§], shown here in Fig. [TJ These two equations 
involve the parameters go (the dimensionless strain rate) , 
Sy (needed in order to convert from dimensionless stresses 
s to measured stresses in units GPa), /x, Co, Xooj and the 
four initial values of the effective temperature x(0) — Xo- 
We can obtain a certain amount of information about 
these parameters directly from features of the stress- 
strain curves before doing any computation. Note from 
Eqs. (|2.7p and (|2.1ip that no plastic deformation occurs 
for s < 1, no matter how large x might be. Thus the 
lowest observed value of a stress that marks a departure 
from elastic behavior is an upper bound for Sy. The bot- 
tom stress-strain curve in Fig. [U which is the curve with 
the largest initial value of ^, seems to have a break point 
that is about a factor of three less than the flow stress s/, 
which we see from the figure is about 1.06 GPa. There- 
fore, Sy ^ S//3 ^ 0.35 GPa, and s/ ^ 3. We take the 
shear modulus = 46 GPa directly from the slope of the 
stress-strain curves in the elastic region; thus we choose 
/2 = 130. Setting the right-hand side of Eq. (|2.9p to zero, 
we find that the flow stress s/ satisfies 

go = 7To = 2eoe-i/x-g(S/). (3.1) 

Using this relation with the DA estimate of 7 and the 
value of Xoo given below, we confirm that tq is of order 
femtoseconds as implied by the Stillinger- Weber interac- 
tions used in these simulations. 

Once we have inserted Eq. (|3.ip into Eqs. (|2.9p and 
(|2.10p , we can solve these equations and vary the remain- 
ing parameters to fit the stress-strain curves in Fig. 15 
(a) of i. We find that the fits are disappointingly insen- 
sitive to our choice of eo exp(— 1/xoo) so long as we stay 
within our constraint that this product be small. How- 
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ever, our later analysis of the liquidlike fraction (j) reveals 
that the internally consistent value of Xoo is 0.065, which 
is slightly bigger than our a priori order-of-magnitude 
estimate based on the experimental data, as shown in Q 
for example, but seems well within the accuracy of our 
approximation for Xoc in [s] and the uncertainties of these 
small-scale simulations. Therefore we have chosen eg = 1 
and Xoo = 0.065. Agreement between theory and the nu- 
merical simulations then can be obtained for all four of 
the stress-strain curves by setting cg = 0.18 and adjust- 
ing only xo for each curve. Our best-fit results are shown 
in Fig. [TJ The corresponding values of xo (from top to 
bottom in the figure) are 0.0559, 0.0580, 0.0612, 0.0680. 
Note that there is a gap in xo between the lowest three 
values, for which the stress-strain curves are peaked, and 
the largest, which shows no peak. 

If it were useful to do so, we could improve the agree- 
ment between the simulated stress-strain curves and our 
theory by making small adjustments of Sy and /i for each 
curve, consistent with the likelihood that the four rela- 
tively small computational systems are not exactly com- 
parable to each other in their as-quenched states. In 
the spirit of the discussion to follow, however, we have 
chosen to assume that the four systems reach effectively 
identical steady states after persistent shear deformation, 
and to attribute the small discrepancies to the statisti- 
cal uncertainties visible in the data. The only systematic 
feature of the stress-strain data that is not recovered by 
the theory is the slow decrease of the flow stress with 
increasing strain in the top curve (the one with the most 
pronounced stress peak). We believe that this behavior 
is caused by the emergence of a nascent shear band as 
illustrated in Fig. 16(b) of Q. 



IV. SECOND THOUGHTS 

Fitting the stress-strain curves, however, is jinly a part 
of the challenge of interpreting the DA data [3] . We also 
must understand how the liquidlike fraction and the 
density p relate to our STZ variables, especially the ef- 
fective temperature x- Therefore, it is essential to know 
whether we can trust the values of xo deduced from the 
STZ analysis of the stress-strain curves. We argue in the 
following paragraphs that these values are not quantita- 
tively reliable. 

One of the key tenets of plasticity theory is that tran- 
sient peaks observed in the stress-strain curves for well 
annealed samples occur because there is an initial lack 
of plasticity carriers in these systems, and new carriers 
must be generated by deformation before plastic flow can 
begin and the stress can relax. Demkowicz and Argon 
use the correlation between stress peaks and the liquid- 
like fraction to argue that </> is a direct measure of 
the population of plasticity carriers. The STZ theory, as 
presently constituted, predicts stress peaks when - and 
only when - the initial STZ density is small. 

This tenet is not confirmed by the DA simulations. 



Note first the data for as a function of strain 7 shown 
by DA in Fig.l5 of [|, also shown here in Fig [51 There 
are four curves. Two of them have small initial values of 
(/)(7 — 0) = (j)o and exhibit stress peaks. The correspond- 
ing values of 0(7) rise monotonically to the steady-state 
value = 0.46 as expected. Another curve starts with 
00 > 0OO- The corresponding stress-strain curve has no 
peak, and ^(7) decreases to 0oo, again in accord with 
expectations. In one case (the one with 0o = 0.46), how- 
ever, 00 is slightly above 0oo but the stress still shows a 
peak. Moreover, the corresponding density p decreases 
with 7, implying that the denser liquidlike fraction is de- 
creasing during the deformation. 

This kind of behavior appears elsewhere in the DA pa- 
pers [1, H, 0, [1] . It is clear in the constant-volume, MD 
simulations, where one of the stress-strain curves has a 
peak, but the corresponding 0(7) remains nearly con- 
stant, and the pressure increases instead of decreasing as 
it should if the liquidlike fraction were growing. And, as 
we note below, all of the PEM simulations in Q show 
stress peaks, even the ones in which 0o is large and com- 
parable to 0OO. 

It seems to us that the most likely explanation for these 
discrepancies is that the ubiquitous stress peaks are arti- 
facts of the step-strain simulations. It is also quite possi- 
ble, of course, that the STZ theory does not adequately 
account for the way in which stresses respond to rapid 
changes in the strain rate, and that the DA simulations 
are simply out of range of our STZ analysis. Perhaps 
both explanations are in part correct. Nevertheless, we 
prefer the first explanation for the following reasons. 

Demkowicz and Argon performed only MD simula- 
tions, and not PEM, at constant-pressure; but a number 
of features of their PEM results may be relevant here. 
They report that cascades of events were evident in their 
PEM simulations but were hard to detect in MD. Mal- 
oney and Lemaitre d, [13], who used only PEM, found 
that cascades were prevalent and sometimes so large that 
they spanned their systems, which were only two dimen- 
sional but larger in numbers of atoms than those used 
by DA. A related feature of the DA PEM simulations is 
that, even for rapidly quenched samples with large initial 
disorder, the stress-strain curves exhibit marked stress 
peaks followed by strain softening. We suspect that these 
behaviors may be generally characteristic of step-strain 
processes, and correspondingly uncharacteristic of con- 
tinuous strain mechanisms that are more common in the 
laboratory. 

Note that, between each strain increment in PEM, the 
system has a probability of dropping into a low energy 
state from which it cannot escape without the applica- 
tion of a large force. The large energy released when 
such a trapped state is destabilized may trigger cascades 
of smaller events. This trapping mechanism may be espe- 
cially important at the beginning of a shear deformation, 
because then the energy minimization starts with a disor- 
dered, as-quenched system. As a result, the first energy 
drops may be large, and the initial stresses required to 
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set the system into motion may be anomalously high. 
Thus we expect transient stress peaks in PEM, even for 
initially disordered systems; and we expect large stress 
fluctuations even in steady-state conditions. That is ex- 
actly what is seen by Demkowicz and Argon. 

The difference between step strains and continuous 
shear has been demonstrated recently in bubble raft ex- 
periments by Twardos and Dennin [llj . Like Demkowicz 
and Argon, these authors subjected their strictly ather- 
mal system to both steady, slow shear and to discrete 
shear steps followed by relaxation periods long enough 
for most motion to cease. They monitored the stress 
during both procedures, with particular interest in the 
distribution of stress drops associated with irreversible 
plastic rearrangements. One of their most interesting re- 
sults is that the average size of stress events decreases 
with decreasing shear rate for continuous strain, but in- 
creases for step strains. That is, the stress relaxes via 
a larger number of smaller drops for continuous defor- 
mation than for step-strain motion. Their interpretation 
of this result seems to us to be roughly consistent with 
our discussion of PEM simulations in the preceding para- 
graph; but neither they nor we claim a full understanding 
of this phenomenon 

In their MD simulations, Demkowicz and Argon use 
step strains as opposed to a continuous strain rate. Ap- 
parently, allowing the system to relax for a time at 
T = 300 K between strain steps removes the effects of at 
least some of whatever structural irregularities are pro- 
ducing cascades and stress peaks in the PEM simulations. 
It seems likely, however, that some elements of the ather- 
mal PEM-like, step-strain behavior persist in the MD re- 
sults. Our hypothetical low-energy trapping states may 
be less likely to be sampled in step-strain MD, in which 
case the initial transients and steady-state stress fluctu- 
ations might be smoother, as indeed they are. But the 
step-strain MD procedure is not the same as a continuous 
one in which the strain is incremented on the time scale 
at which the molecular motions are resolved. In the DA 
simulations, the strain is incremented only once in about 
ten or more atomic vibration periods. We see little rea- 
son to believe that these two simulation procedures will 
produce precisely the same responses to rapidly changing 
driving conditions. 

This situation forces us to conclude - reluctantly in 
view of the quality of the apparent agreement between 
the simulations and STZ theory shown in Fig. [T]- that 
the values of the xo's stated above and in the caption to 
Fig. [1] cannot be trusted. 



V. QUASITHERMODYNAMIC HYPOTHESES 

Our inability to deduce the xo's from the stress-strain 
data has led us to probe more deeply into the mean- 
ing of the effective temperature in glass dynamics. We 
base our analysis on a set of quasithermodynamic hy- 
potheses in which we assume that the steady-state prop- 
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FIG. 2: Upper panel: The equation of state (p — 0o(x) 
obtained by fitting the four data points (filled circles) plus 
4>oo ~ <^a(xoo) ~ 0.458 (filled square). The 0o's are the 
DA numerical simulation data [3|; the adjusted xo's are 
0.0559, 0.0589, 0.0663, and 0.0749, with = 0.065. The 
solid line is our fit to these points in Eq. (|5.ip . The open 
circles are the observed (j>o's plotted against the unadjusted, 
best-fit, Xo's. The dashed line is intended only as a guide 
for the eye. Lower panel: The equation of state p — po(x) 
obtained by a procedure similar to the one described above 
except that the point poc ~ 1.0335 (filled square) is not in- 
cluded in the fit. We attribute its displacement from the 
equation of state to a stress-induced dilation. The density is 
normalized by pc, the density of crystalline silicon (diamond 
cubic). The solid line is our smooth fit to the data. Again, 
the open circles and and dashed line indicate the unadjusted 
equation of state. 



erties of the configurational degrees of freedom below the 
glass transition are determined by the effective temper- 
ature Teff — (EsTZ /ks) ill close analogy to the way 
in which they are determined by the bath temperature 
T above that transition. In true thermodynamic equilib- 
rium, extensive quantities such as the density p, the inter- 
nal energy U, or the liquidlike fraction (j) in the DA sim- 
ulations are functions of T and the pressure P. In other 
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words, these quantities obey equations of state. Below 
the glass transition, the configurational degrees of free- 
dom - that is, the positions of the atoms in their inherent 
states - fall out of equilibrium with T because thermally 
activated rearrangements are exceedingly slow or impos- 
sible. The most probable configurations in such situa- 
tions maximize an entropy, thus the statistical distribu- 
tion of these configurations is Gibbsian with T = T^ff. 

Accordingly, our first hypothesis is that the high-T 
equilibrium equations of state are preserved in the glassy 
state as expressions for the configurational parts of the 
extensive quantities p, U , (j), etc. in terms of the intensive 
quantities T^ff, P, and the shear stress s. In particular, 
we propose that the effective temperatures Xoi multiplied 
by EsTZ /ks, are the bath temperatures Tq at which the 
DA simulation samples fell out of thermal equilibrium 
during cooling, and that the configurational parts of p, 
U, (j), etc. were fixed at those temperatures. Our special- 
ization to configurational parts recognizes that, for exam- 
ple, p undergoes ordinary thermal expansion at small T 
and that the total internal energy U includes the kinetic 
energy. Both of those non-configurational parts are un- 
interesting for present purposes. Since we consider only 
P = situations and work only at fixed T, we omit ex- 
plicit dependences on those variables. Shear dilation may 
be relevant, especially for determining p; but so long as 
we are considering only pre-deformation properties with 
s = 0, we also omit explicit s dependence. We return 
later to the question of shear dilation, and we also post- 
pone a discussion of the potential energy. For the mo- 
ment, we write simply: 4> — 0o(x): P — Poix)- 

As a first test of this hypothesis, we show by the open 
circles and dashed lines in both panels of Fig. O the 
functions 0o(x) a-i^d po(x) obtained with the measured 
values of and po and our best-fit values of xo from Fig. 
[TJ We do, in fact, find qualitatively plausible behavior. 
With our uncertainties about the Xo's, however, we are 
obliged to look harder and bring other considerations to 
bear on this analysis. 

We know some qualitative features of these equations 
of state. First, we are dealing with a glass transition, 
and therefore we expect that 4>o{x) should extrapolate 
to some fixed value, possibly zero, at a nonzero Kauz- 
mann temperature. In the following discussion, we de- 
fine xk to be the effective temperature at which the ex- 
trapolated 00 (x) vanishes; but we make no assumption 
about whether this is actually a Kauzmann temperature, 
or whether an ideal glass transition actually occurs at this 
or any other point. 

Another feature is suggested by the binary, solid- 
like/liquidlike, nature of the atomic structure described 
by Demkowicz and Argon. If this were just a simple 
binary mixture, and if the energy of the mixture were 
proportional just to the density of the higher energy, liq- 
uidlike component, then the entropy would be a maxi- 
mum at (f) = 1/2 and the temperature would diverge at 
that point. Obviously, amorphous silicon is not so sim- 
ple. Nevertheless, the observed (f> does seem to have an 
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FIG. 3: Modified stress-strain curves using xo = 
0.0559, 0.0589, 0.0663, and 0.0749. All other parameters are 
the same as those used in Fig. [1] 

upper bound not too far above — 1/2; thus, we expect 
that 00 (x) saturates at some value 0o ~ 0iim > 0.55 for 
large x- 

Yet another consideration is that, in principle, we 
ought to be able to deduce the To's directly from the 
graphs of density versus temperature for different quench 
rates shown in Fig. 1 of 3]. The slowest quench is the 
closest to full equilibrium, and the points at which the 
faster quenches fall away from it are the three higher 
Tq's. The lowest To must be near the inflection point on 
the slow-quench curve. Unfortunately, the best we can 
do with the data at hand is estimate the To's to about 
±20 K, and even then we must be cautious about sys- 
tematic uncertainties in the simulations. Nevertheless, 
we can use this process as a rough check on the modifi- 
cations of the xo's that we propose below. 

We come now to our second and more speculative 
hypothesis - that the point (xoo, 0oo) ought to lie on 
the equilibrium equation-of-state curve; that is, (poc = 
0o(Xoo)- Here we are assuming that persistent shear de- 
formation in an athermal system rearranges the atomic 
configurations in a way that is statistically equivalent to 
thermally driven rearrangement, except that the relevant 
temperature is the effective disorder temperature Tg/ / in- 
stead of the bath temperature T. We also are assuming 
that 0OO, being a ratio of two populations, is insensitive 
to changes in the volume of the system that might occur 
during constant-zero-pressure simulations. In contrast. 
Poo would decrease if, for example, the system undergoes 
a shear-induced dilation. 

With these considerations in mind, our next step is 
to refine the estimates of the To's by requiring that 
they produce smooth equations of state with the qualita- 
tive features hypothesized above, and that the observed 
(j)oo = 0.458 be equal to 0o(Xoo)- Accordingly, we have 
computed 4> = 4>o{x) p = po(x) by adjusting the 
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Xo's so as to optimize agreement with all the previously 
discussed constraints. Our results are shown by the solid 
lines in Fig. [21 As expected, poo < Po(Xoo), because the 
flow stress s/ at that point seems easily large enough to 
cause dilation. (Throughout this paper, values of p are 
in units of pc = 2323.8 kg/m , the density of crystalline, 
diamond-cubic, silicon.) The effective Kauzmann tem- 
perature is Xk = 0.051, and the upper bound for is 
(f'lim = '/'(x^c») — 0.6. The smooth curve that we have 
fit to the equation of state for 0, shown in the upper 
panel of Fig. [2l is 

0o(x)=0iim (l-e-'^f™)), (5.1) 

where a = 100. Other important parameters are the 
four To's: 1100 ivT, 1160 i^, 1305 /C, 1475 which are 
consistent with our crude direct estimates from Fig. 
1 in [3]. Note that the large gap between the low- 
est three unadjusted To's and the upper one has disap- 
peared, and that the adjusted points fit on a smoother 
curve than the unadjusted ones. We also find that 
Too = 1280 iC. Then, with Xoo = 0.065, we obtain 
EsTz/ks — 21,000-?^, or about 1.3 ev. The associated 
values of xo are 0.0559, 0.0589, 0.0663, and 0.0749. 

Using these adjusted values of the Xo's, we have re- 
computed and plotted the stress-strain curves in Fig. [3l 
shown again in comparison with the simulation data. 
As expected, the non-STZ peak has disappeared, and 
there is more rounding in the bottom, most rapidly 
quenched case, but the stress peaks for the two most 
deeply quenched systems remain almost unchanged. 

Finally, within the context of quasithermodynamic 
equilibrium, we note that a relation between the poten- 
tial energy and x has been proposed and confirmed by 
Shi et aL.[l2i| These authors describe molecular-dynamics 
simulations of shear banding in two-dimensional, non- 
crystalline, Lennard- Jones mixtures. Their analysis of 
these simulations goes beyond our own in at least one 
important way. They note that the plastic strain rate 
predicted by the STZ theory, as well as by other fiow- 
defect theories, has the form exp(— 1/x) multiplied by a 
function of the shear stress, essentially our factor q{s) 
defined in Ea. (|2.1ip . Force balance requires that s be a 
constant across the shear band, thus a measurement of 
the position-dependent strain rate is a measure of the 
position dependence of x- Shi et al. then postulate 
that the potential energy depends linearly on x- They 
compute the position-dependent potential energy directly 
from their simulation data and find that it maps accu- 
rately onto X as predicted by the spatially varying strain 
rate. The importance of this observation is that they are 
looking at a weakly nonequilibrium situation in which the 
potential energy and x are varying continuously in space, 
but only very slowly in time, along the postulated linear 
equation of state. Their system remains in quasithermo- 
dynamic equilibrium during this variation because, unlike 
the DA simulations discussed here, it is not undergoing a 
fast transient response to a sudden change in the applied 



stress. Thus their result anticipates and fits accurately 
into our quasithermodynamic picture. 

VI. DEPARTURES FROM 
QUASITHERMODYNAMIC EQUILIBRIUM 

The next question is whether the equilibrium equations 
of state shown in Fig. [51 are obeyed under nonequilibrium 
conditions during deformation. Apparently they are not. 
To sec this, in Fig. [H we have added to our equation-of- 
state graphs four data sets showing the observed values 
of (j) and p as functions of our calculated values of x for 
the four different deformation histories. In all cases, the 
observed values of (p and p initially fall off the equilib- 
rium curves and do not come together again until they 
approach Xoo ■ This nonequilibrium behavior can be seen 
more directly by comparing the three panels of Fig. 15 
in [3J. It is clear there that the rate at which the stress 
relaxes to the flow stress is faster than the rates at which 
either cf) or p relax to their steady-state values. 

To account for the transient, nonequilibrium behavior 
of (j), we propose an equation of motion analogous to Eq. 
(17) in [5|, but with a form similar to our Eq. (|2.10p : 

^ = ^e-'/^rsqis){Mx)-<l^), (6.1) 

d'j ci qo 

where ci is a constant similar to Cq. The idea here is that 
the disorder described by x may rise at its own rate, but 
the liquidlike-solidlike reorganization associated with 
may not catch up instantaneously. We assume that the 
underlying mechanism that determines these rates is still 
the rate of energy dissipation. There is no other compara- 
bly simple and basic coupling between mechanical defor- 
mation and the internal degrees of freedom that satisfies 
the requirement that it be a non-negative scalar. Note 
that Ea. (|6.ip is supplementary to the STZ equations of 
motion, Eas. (|2.9[ ) and (|2.10p : it assumes that the inten- 
sive quantities 5(7) and xil) arc unchanged from the STZ 
predictions and that they continue to control the behav- 
ior of (j) even away from equilibrium. This is a strong 
assumption, especially near the initial stress transients 
where we already know that there is some mismatch be- 
tween the STZ theory and the DA simulations. 

Our results for the four functions 4>{x)i determined 
using Eq. (|6.ip . are shown in Fig[5l The corresponding 
functions 4'{'y) are shown in Fig[Sl here in comparison 
with the DA data. The only adjustable parameter is ci, 
which we choose to be 0.42. The agreement seems to be 
within the uncertainties of the data. 

We continue our development of the nonequilibrium 
quasithermodynamic theory by writing an equation anal- 
ogous to Eq. (|6.ip for the function p(7): 

^ = i£l e-Vx s q(^s) (Mx) -P- Spo) , (6.2) 

Here, dpo = 0.012 is the dilation-induced shift of the equi- 
librium density shown by the displacement of po^ from 
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FIG. 4: Upper panel: The DA numerical simulation data for 
versus the theoretical x for the four cases considered in Ref. 
[31 . The equation of state = 4>o{x) '^^^ added to stress the 
nonequilibrium nature of the deformation-induced dynamics 
of 0. Lower panel: The corresponding plot for p/pc- The 
points (Xoo, <i>oo) and (xoo, Poo) are marked by solid squares. 



the equilibrium curve po (x) i^i the lower panel of Figl^l 
To integrate Eq. (|6.2p . we have approximated po(x) by 
a smooth polynomial. The factor ci — 0.42 needed to 
fit the density data is the same as in Eq. (|6.ip . The 
corresponding graphs of yo(7), along with the DA data, 
are shown in Fig. [T] Again, the agreement seems to 
be within the uncertainties; but here we are pushing the 
theory too far for comfort. Our quasithermodynamic hy- 
potheses imply that po(x) ought to be a function of the 
shear stress s as well as x, and that the dilation approx- 
imated here by 5pQ should be part of that generalized 
equation of state. In particular, 5pQ should be the dila- 
tional change in the density when s is equal to the flow 
stress s / . We have in fact tried 5pQ oc (as in nonlinear 
elasticity), but the results are distinctly unsatisfactory at 
small 7 where neither Eq. (|6.2p nor the STZ theory itself 
may accurately describe the fast transient. Equations 
(16. ip and (|6.2p do account for the relatively slow relax- 
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FIG. 5: Upper pa,ne\:<j) as a function of x using the nonequi- 
librium Eq. (|6.ip with ci = 0.42. The equilibrium equation of 
state was added for illustration. The curves should be com- 
pared to the DA data shown in the upper panel of Fig. |4] 
Lower panel: p as a function of x using the nonequilibrium 
Eq. (16. 2p with ci = 0.42. These curves should be compared 
to the DA data shown in the lower panel of Fig. |4] The points 
(Xoo, 0oo) and (xoo, Poo) are marked by solid squares. 



ation of 0(7) and p(7) as compared to that of 3(7). This 
level of success seems to be as much as we can expect 
from the theory at this stage in its development. 



VII. SUMMARY AND CONCLUDING 
REMARKS 

The unique aspect of the Demkowicz and Argon simu- 
lations is their measurement of extensive quantities - the 
mass density p and the liquidlike fraction (/) - in parallel 
with conventional stress-strain curves. They argue con- 
vincingly and importantly that (j) is closely related to the 
density of plasticity carriers. However, their remarkably 
complete report includes some results that make it seem 
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FIG. 6: Theoretical predictions of 0(7) for all four quenches, 
based on Eq. (|6.ip with ci = 0.42, compared with the DA 
simulation data. 
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FIG. 7: Theoretical predictions of ^(7) for all four quenches, 
based on Eq. (16. 2|) with ci = 0.42, compared with the DA 
simulation data. 



that this relationship may be neither simple nor direct. 
The main purpose of our investigation has been to learn 
more about that relationship. 

Our proposed interpretation of the DA simulations is 
based on the athermal STZ theory [1], which we believe 
captures the central features of amorphous plasticity - 
at least for processes that are not too rapidly varying as 
functions of space or time. In order to discuss quanti- 
ties like (j), however, we have had to go beyond the STZ 
theory. We have hypothesized that, in low-temperature, 
steady-state, nonequilibrium conditions, p and (f) and pre- 
sumably other such quantities are related to the effective 
disorder temperature x ^^o, quasithermodynamic equa- 
tions of state. With this hypothesis, we have developed 



a theoretical interpretation of the DA simulations that 
seems physically satisfying, internally self-consistent, but 
interestingly incomplete. If confirmed by further tests, 
this quasithermodynamic theory could become a useful 
tool for predicting the nonequilibrium mechanical behav- 
ior of amorphous solids. 

There are many open issues. Perhaps the most ur- 
gent of these is the interpretation of the DA step-strain 
simulation technique. The DA data do show that a well 
annealed sample with an initially small <f) exhibits a tran- 
sient peak in its stress-strain curve. However, the con- 
verse seems not necessarily to be true. Stress peaks some- 
times appear in the DA results in cases where the initial 
(j) is large and does not increase during plastic deforma- 
tion, a behavior that is inconsistent with the presumed 
relation between and the density of plasticity carriers. 
After examining other possibilities and looking at other 
examples of step-strain procedures (see Sec lIV[) . we have 
concluded that the stress peaks observed in cases with 
large initial (j) &re most likely artifacts of the numerical 
step-strain procedure. Once we allowed ourselves that 
flexibility and found alternative ways of estimating the 
initial effective temperatures for those cases, we found 
that our quasithermodynamic picture fits together quite 
well. 

In our opinion, one of the most important next steps 
in determining the limits of validity of the quasithermo- 
dynamic STZ theory would be to redo the DA analysis, 
including the measurement of the liquidlike fraction ^, 
with continuous strain MD. So far as we know, nobody 
before Demkowicz and Argon has ever done anything like 
this " making independent, simultaneous measurements 
of the density of plasticity carriers and the stress-strain 
response for differently quenched systems. Discovering 
the relation between continuous and step-strain simula- 
tions in this context seems likely to be extremely useful 
at least for understanding the numerical simulations. It 
could also be a big step forward in understanding amor- 
phous plasticity. 

Another obviously open issue is the validity of the STZ 
theory in situations where the system is responding to 
rapid changes in loading conditions, as in the DA numeri- 
cal simulations where the strain rate is large and is turned 
on instantaneously both at the beginning of the process 
and at each strain step. The present version of STZ the- 
ory |l| is a mean-field approximation in which interac- 
tions between zones are included only on average, and 
fluctuations are neglected. It is possible that neither the 
STZ theory as presently formulated, nor the conventional 
explanation of stress peaks in terms of plasticity carri- 
ers, can fully account for stress transients during rapid 
changes in loading conditions, even in continuous-strain 
MD simulations. From this point of view, a continuous- 
strain version of the DA simulations seems doubly impor- 
tant. If new simulations showed no qualitative discrep- 
ancy between continuous and step-strain procedures, we 
would know that the STZ theory is missing some essential 
ingredients. We might then try to include correlations 
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and fluctuations by developing a more detailed statisti- 
cal theory of STZ's with varying thresholds, perhaps with 
noisy interactions between them as proposed recently by 
Lemaitre and Caroli fl^. 

In a similar vein, wc must ask about the limits of valid- 
ity of our quasithcrmodynamic hypotheses. It will be es- 
pecially important to look harder at Eqs. (|6.ip and (|6.2p . 
which determine the rates at which and p relax to their 
quasithcrmodynamic equilibrium values. Our quasithcr- 
modynamic picture is based on the assumption that the 
effective temperature x is the intensive variable - the 
thermodynamic force - that drives the extensive quanti- 
ties 0, p, and the STZ variables A and m. According to 
our equations of motion, A and m are tightly slaved to x- 
We have postulated equations of state relating cj) and p to 
X under equilibrium or steady-state conditions, and have 
further proposed in Eqs. (|6.ip and (|6.2[) . on what seem to 
us to be general grounds, that 4> a-nd p are more loosely 
slaved to x than are A or to during excursions from equi- 
librium. Here, as in the questions regarding the stress 
response, we need to learn whether the equations of mo- 
tion for (j) and p are accurate when those excursions from 
equilibrium are faster, say, than the relaxation rate of x- 

In short, we are asking how far this theory can be 
pushed. Can it, for example, always be used to predict 
stress-strain transients under realistic experimental con- 



ditions, where strain rates are very much smaller than 
those in MD simulations? Or does it generally become 
inaccurate near stress peaks? Can it be used to pre- 
dict plastic deformation near the tip of an advancing 
crack? More generally: Is the effective temperature x re- 
ally such a dominant state variable? What other internal 
variables might become relevant for describing fast pro- 
cesses? What real or computational experiments might 
help to answer such questions? 
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